Pushing the hyperpolarizability to the limit 
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We use numerical optimization to find a one-dimensional potential energy function that yields 
the largest hyperpolarizability, which we find is within 30% of the fundamental limit. Our results 
reveal insights into the character of the potential energy functions and wavefunctions that lead to 
the largest hyperpolarizability. We suggest that donor-acceptor molecules with a conjugated bridge 
with many sites of reduced conjugation to impart conjugation modulation may be the best paradigm 
for making materials with huge hyperpolarizabilities that approach the fundamental limit. 
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The study of nonlinear optical materials has 
been driven by the quest for making the largest 
possible nonlinear response in a medium for use 
in a broad range of optical applications such as 
telecommunications, optical data storage, three- 
dimensional nano-photolithography,0 Q and making 
new materials for novel cancer therapies. Q 

A natural questions arises whether there is a limit to 
the nonlinear response and if so, how does one make a 
molecule to get to the limit. The former has been an- 
swered in the affirmative using sum rules. 1^. 1^. fToL ITU 
IT^ I and the fundamental limit is given by, 
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where N is the number of electrons and E w the energy 
to the first excited state. This result has not yielded 
clear guidance on how to design better molecules, but 
rather has set the bar for assessing the performance of 
molecules. In fact, the largest nonlinear susceptibilities 
of the best molecules fall short of the fundamental limit 
by a factor of 10 3 / 2 .[H[H[3 A Sum-Over-States (SOS) 
calculation of the hyperpolarizability^^ using the ana- 
lytical wavefunctions of the clipped harmonic oscillator 
yields a value that is within a factor of 2 of the limit. |14| 
Thus, the factor-of-thirty gap between the fundamental 
limit and the apparent upper bound of the best molecules 
is not of a fundamental nature. So, it should be possible 
in principle to make materials with second-order suscep- 
tibilities that are a factor of 30 bigger than the currently 
best materials. In this letter, we report on the charac- 
ter of the potential energy function of a one-dimensional 
system that yields a hyperpolarizability, /3, that is within 
30% of the fundamental limit, and use this potential to 
propose a new paradigm for molecular engineering. 

There are two equivalent expressions for p. The stan- 
dard one, flsos, as calculated by Orr and Ward.p^f and 
the dipole-free expression, (3df-^\ The standard one 
is overspecified in the sense that it is possible to pick 
unphysical values of the energies and matrix elements, 
which violate the precepts of quantum mechanics. The 
dipolc free expression, in contrast, is simplified into a re- 
duced form that contains no dipole terms. Since the two 
expressions should yield the same results if all the matrix 
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FIG. 1: Evolution of (a) /3 and too; and (b) function 
of number of iterations of optimization. 



elements are consistent with the sum rules, deviations be- 
tween them can be used as a convergence test.0] How- 
ever, we later discuss how to treat exceptions where (3df 
is accurate and (3sos is not. 

Our approach is to numerically calculate j3 for a spe- 
cific potential, then, to use an optimization algorithm 
that continuously varies the potential in a way that max- 
imizes p - using the convergence test to determine if 
the result is reliable. Our code is written in MATLAB. 
For each trial potential we use a quadratic finite element 
method |17| to approximate the Schrodinger eigenvalue 
problem and the implicitly restarted Arnoldi method 0] 
to compute the wave functions and energy levels. To op- 
timize P we use the Nelder-Mead simplex algorithm [ljj . 

Figure shows the evolution of the hyperpolariz- 
ability as a function of the number of iterations of the 
optimization algorithm, applied to Pdf- All hyperpo- 
larizabilities are normalized to the fundamental limit. 
The initial potential energy function is a hyperbolic tan- 
gent with infinite potential at the boundaries. We select 
this potential because it is localized near the origin, yet 
flat elsewhere, allowing the optimization process to work 
without any initial biases. Also, it meets the criteria of 



2 



being noncentrosymmetric, as required for nonvanishing 
hyperpolarizability. After 7,000 iterations, the algorithm 
converges to just over (3df/(3max = 0.72, the largest 
hyperpolarizability seen to date. No other starting po- 
tential, including polynomials, and fractional exponents, 
leads to larger (3df- 

As a second convergence test, we use the fractional 
deviation from the ground state sum rule, which in terms 
of kqo a s defined in the literature, is given bv.|20| 
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where E n o is the energy difference between state n and 
0, x nm is the position matrix element between state n 
and m, and x r {^ ax is the fundamental limit of x\o- Too ap- 
proaches zero when the optimization process yields max- 
imum Pdf, supporting the fact that the numerical com- 
putations accurately represent the system. When, on the 
other hand, (3sos is optimized, Too is large and the two 
forms of (3 diverge appreciably, illustrating the robustness 
of the dipole-free form and the pitfalls of using flsos for 
optimization. Our calculations of (3 include a total of 15 
states, so a large value of r o is more of an indication 
that not enough states are being used in flsos rather 
than inaccuracies of the wave functions. The fact that 
(3 sos and (3df do not converge illustrates the need for 
two independent convergence tests. 

Figure \\]p shows the evolution of the position matrix 
elements. When (3df is small, many of these matrix ele- 
ments are non-negligible, so many states will contribute 
to (3. By a thousand iterations when the hyperpolariz- 
ability approaches its largest value, all matrix elements 
vanish except for two of them. This clearly shows that 
as the fundamental limit of [3 is approached, the system 
collapses to a three-level model for (3. This is consistent 
with the three-level ansatz, previously proposed, that a 
three-level model describes a system with a hyperpolar- 
izability at the fundamental limit. j^J There is mount- 
ing evidence that the three-level ansatz also is obeyed in 
planar molecules. [23] We stress that our results do not 
prove the converse; that all molecules with hyperpolar- 
izabilities at the fundamental limit must be represented 
by a three-level model. However, since (1) the three- 
level ansatz is used to determine the fundamental limits, 
and, all reliable measurements and calculations yield hy- 
perpolarizabilities below the fundamental limit; and (2) 
numerical optimization yields a three-level model, this 
form of induction provides strong support for the ansatz. 

Figure [21 shows the evolution of the potential en- 
ergy function for a length scale on the order of a large 
molecule. Optimization of the hyperpolarizability clearly 
favors a potential with large oscillations. Interestingly, all 
starting potentials we have studied, independent of their 
initial form - such as polynomials, power laws, and piece- 
wise continuous functions, develop such wild oscillations 
of about the same period when the hyperpolarizability is 
maximized. We find that these oscillations serve to local- 
ize each eigenfunction at different positions. Early in the 
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FIG. 2: Evolution of potential energy as a function iteration. 
Inset shows a class of proposed molecules, where refers 
to a continuation of the bridge using the same theme. 



- Ground State Wavefunction 



Excited State Wavefunction, n 




10 12 14 16 18 20 
Position 



FIG. 3: Potential energy and energy eigenfunctions after 
7,000 iterations, when /3df is optimized. 



optimization process, the wavefunctions are all delocal- 
ized while at 7,000 iterations, most of the wavefunctions 
are mutually non-overlapping. 

Figurc[3]shows the eigenfunctions and potential energy 
function after 7,000 iterations. Only excited states num- 
ber 1 and 6 have appreciable overlap with the ground 
state and with each other. As such, the only term that 
contributes to the hyperpolarizability is proportional to 
XmXiQX§Q / E\qEqq. Note that with other starting poten- 
tials, states other then 0, 1, 6 may be important. Most 
significantly, (3 is maximum for a three-level system. 

Under the three-level ansatz, the normalized hyperpo- 
larizability can be expressed as a product of two func- 
tions, p/0MAxf(E)G(X)M where X = x w /x% AX 
and E = E 10 /E 2 q. "1" and "2" label the states with 
the largest transition moments to the ground state. The 
function f(E) is maximum when E = (/(0) = 1) and 
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G{X) is maximum when X = (G( "n/3) = 1). So, 
a hyperpolarizability that approaches the fundamental 
limit should have two dominant states that have well 
separated eigenenergies; and, xio ~\fixiQ AX . The op- 
timized wavefunctions yield £'0.314 and X — 0.775, 
so f(E) = 0.892 and G(X) = 0.999. The transition 
moment is near optimal but the energy level spacing 
can be improved. If the system were truly 3 levels, 
P/0 MAX = f(E)G(X) = 0.89 compared with the op- 
timized value of 0df I Pmax = 0.72. 

One might argue that our results are artificial because 
we are using only 15 states, and, because our space is 
bounded by an infinite potential. These two issues go 
hand in hand. In our calculations, the potential energy 
function is represented by a cubic spline of 20 points. 
As such, it would be numerically impossible to localize 
more than about 15 eigenfunctions. Furthermore, higher- 
energy wavefunctions would interact with the walls; and, 
for high-enough energies, the wavefunctions would have 
the character of a particle in a box. So, to increase 
the number of energy levels, one would need to increase 
the number of points in the spline, which would lead to 
more oscillations that would localize more of the excited 
states. Based on the pattern that we have observed, 
when increasing the number of states and making the 
potential correspondingly broader, we expect that more 
non-overlapping wavefunctions would develop, yielding 
a similar conclusion. So, in the limit of infinitely wide 
space, the conclusion should be the same. 

In summary, we have applied an optimization process 
to determine what potential energy function yields the 
largest value of the hyperpolarizability. As a function of 
iteration, the hyperpolarizability approaches the funda- 
mental limit, but never exceeds it. When (3 converges 
to the maximal value, the expression for (3 collapses to 
a three-level model, supporting the three-level ansatz. 



Strong oscillations in the potential energy function serve 
to isolate the wavefunctions to prevent overlap between 
all states but three of them. This behavior suggests that 
long linear organic molecules with regions of conjugation 
(representing the dips in the potential) interspersed by 
regions of reduced conjugation (representing the peaks) 
and flanked by a donor and acceptor to break the sym- 
metry, would be a promising new paradigm. 

The inset in Figure [21 shows one example of a new 
molecular paradigm that should exhibit the required 
modulation of conjugation over a single and double bond 
pair spatial period. The dashed bonds represent the ab- 
sence or presence of rings to control the degree of con- 
jugation and/or molecular planarity. Note that at 100 
iterations, the oscillations are on the order of leV (com- 
mon energy scale for molecules), where 0/ 0maxO-2 ~ well 
above the apparent limit of 0/0max = 0.03. So, our 
new paradigm should lead to a tenfold improvement over 
today's best molecules. Alternatively, our method could 
guide the fabrication of multiple quantum well structures 
using a variety of layered organic molecules that impart 
to it an oscillating potential energy profile. 

We plan to apply our general method to higher-order 
nonlinear-optical susceptibilities. In addition, it can be 
applied to resonant processes and to any general com- 
binations of laser wavelengths. With regards to the sec- 
ond hyperpolarizability, 7, since past calculations suggest 
that the three-level ansatz also should hold, we would ex- 
pect the same sort of conclusion. Should this approach 
be successful, new and exceptional materials could be 
developed to make more efficient devices and novel ap- 
plications. 
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